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I.  INTRODUCTION 


A  line  by  line  calculation  is  the  most  accutate  and  general  means  for 
computing  the  atmospheric  transmission  of  radiant  energy.  With  the  aid  of  a 
modem  computer  and  the  recent  convenient  availability  of  tables  of  line  para¬ 
meters  stored  on  magnetic  tape,^  such  calculations  can  now  be  routinely 
carried  out. 

For  many  applications,  the  quantity  of  interest  is  the  mean  value  of  the 
transmittance  over  a  band  Interval  Av 

T  -  /,  T(v)  dv  (1) 

Av  •'  Av 

The  total  numerical  effort  required  to  compute  T  using  the  line  by  line  method 
Is  quite  large  and  any  extensive  calculation  of  mean  transmittance  values  can 
be  quite  time  consuming,  even  on  a  fast  computer. 

Band  model  theory,  on  the  other  hand,  provides  a  simple  parameterized 
formula  for  computing  the  mean  band  transmittance,  which  is  much  more  effi¬ 
cient  than  the  line  by  line  calculations.  However,  since  it  is  based  on 
certain  simplifying  assumptions  It  must  be  tested  for  accuracy.  Previous 
studies  comparing  the  statistical  band  model  with  line  by  line  calculations 
have  shown  a  reasonable  agreement  between  the  two  methods.  This  comparison, 
however,  was  carried  out  for  small  to  moderate  optical  path  distances  where 
the  transmlttances  were  In  the  range  T  *  >  0.1.  Comparisons  have  not  been 
made  In  the  very  long  path  length  regime  where  f  <  0.1.  In  this  report,  we 
present  the  results  of  such  a  comparison.  The  details  are  given  In 
Section  II. 

We  will  conclude  from  this  comparison  that  the  band  model  does  not.  In 
general,  give  very  good  results  in  the  regime  where  T  ■  <  0.1.  In  fact,  as 


^L.  S.  Rothmon,  Appl.  Opt.  20 ,  791  (1981). 

^A.  Goldman  and  T.  Kyle,  Appl.  Opt.  _7_,  1167  (1968). 


will  be  seen,  the  transmittance  can  be  In  error  by  several  orders  of  magnitude 
and  Is  shown  to  have  the  wrong  asymptotic  behavior  as  the  optical  path  length 
approaches  Infinity. 

In  Section  III,  we  derive  new  parameterized  formulas  for  computing  the 
mean  band  transmittance  T(x).  These  formulas  are  accurate  for  all  values  of 
optical  distance  x.  Including  the  limit  as  x  approaches  Infinity.  The  results 
of  comparisons  with  both  the  precise  line  by  line  calculations  and  statistical 
band  model  calculations  will  also  be  presented. 

Section  IV  Is  a  summary  and  discussion. 


II.  COMPARISON  OF  STATISTICAL  BAND  MODEL  AND  LINE  BY  LINE  TRANSMITTANCE 


The  line  parameter  data  for  these  calculations  were  obtained  from  the  Air 
Force  Geophysics  Laboratory  (AFGL)  line  atlas. ^  This  compilation  provides 
data  for  seven  gas  species:  H2O,  CO2,  O3,  N2O,  CH^,  CO,  and  O2.  The  data 
given  for  each  line  include  the  line  position  (cm“^),  the  line  strength  S^ 
(cm  ^/molecule  -  cm~^  at  296‘’K),  the  line  half-width  (cm~^  for  air  broaden¬ 
ing  at  296®K  and  1  atm  pressure),  and  the  energy  of  the  lower  level  of  the 
transition  E^  (cm~^).  The  subscript  i  labels  the  line,  while  the  superscript 
0  on  the  strength  S^  and  the  half-width  indicate  that  these  quantities  are 
taken  directly  from  the  atlas  and  are,  therefore,  appropriate  for  a  pressure 
of  1  atm  and  temperature  of  296“K. 


The  line  strength  and  the  width  parameter  y^  at  any  other  temperature 
and  pressure  are  computed  from  and  y^  using  the  same  formulas  employed  by 
C.  M.  Randall  in  his  general  line  by  line  computer  program,  INHOM.  A  discus¬ 
sion  of  these  formulas  is  given  in  Ref.  3. 


All  calculations  have  been  carried  out  assuming  a  Lorentz  pressure  broad¬ 
ened  line  shape  function.  The  spectral  absorption  coefficient  is  given  by 


k(v) 


"i  h  ^i 


(v  -  -  y^^ 


(2) 


The  quantity  n^^  is  the  density  of  the  gas  for  which  and  y^  characterize  a 
line.  The  units  of  n^  used  in  this  study  are  molecules /(cm^  -  km),  i.e.,  the 
number  of  molecules  in  a  column  1  km  long  and  1  cm^  in  cross  sectional  area. 
By  including  the  density  factor  n^^  in  the  definition,  k(v)  has  the  convenient 
dimension  1/km.  The  transmittance  at  wave  number  v  is  given  by 


T(v,  x)  *•  exp  l-k(v)  xj 


(3) 


3 

S.  J.  Young,  Band  Model  Parameters  for  the  2.7-pm  Bands  of  H2O  and  CO^  in  the 
100  to  3000“K  Temperature  Range,  TR-0076(6970)-4 ,  The  Aerospace  Corp. 

(31  July  1975). 
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where  x  Is  the  distance  measured  In  kilometers.  The  mean  band  transmittance 
Is  then  computed  by  averaging  T(v,x)  over  the  band  Interval  Av  [see  Eq.  (1)] 


T(x)  -  —  /  exp  l-k(v)  x]dv 

Av  \ 

Av 


(4) 


In  the  statistical  band  model  with  exponential-tailed  Inverse  line 

4  3  4 

strength  distribution,  the  mean  band  transmittance  Is  given  by  ’ 


T(x) 


exp 


-M] 


(5) 


where  k,  and  y  three  band  model  parameters.  The  band  model  para¬ 
meters  can  be  expressed  as  certain  averages  of  the  line  parameters  over  the 
spectral  Interval  Av.  The  parameter  k  Is  ieflned  by  the  expression 


k 


1  ^ 
IT  ^ 

AV  i.j 


(6) 


where  L  Is  the  number  of  lines  In  the  Interval  Av*  The  parameter  y  Is  the 
average  line  width 


Y 


Yl 


(7) 


and  6^  Is  a  measure  of  the  effective  average  distance  between  lines  In  Av 


^  “  r^-  "l  ^1  ^1^ 


k  y  1-1 


(8) 


A.  HqO  line  parameters 

The  model  system  used  for  the  calculations  In  this  report  has  the  follow¬ 
ing  specifications:  pressure,  1013  mbar;  temperature,  300‘’K;  concentration  of 
H2O,  0.026;  concentration  of  O2,  0.21.  The  wave  number  range  In  which  the 
calculations  are  carried  out  Is  2800  to  3400  cm  (Only  the  H2O  molecules 


^W.  Malkmus,  J.  Opt.  Soc.  Am.  57 ,  323  (1967). 
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have  absorption  lines  in  this  range;  the  ©2  concentration  is  specified  because 
of  its  influence  on  the  width  parameters  of  the  H2O  lines.)  The  density  n  of 
H2O  molecules  in  this  model  system  is 

n  "  6.36  X  10^^  i^olecule£ 
cm  -  km 

The  spectral  band  Interval  is  chosen  to  be  Av  ■  20  cm~^.  A  band  will  be 
identified  by  its  mean  wave  number.  Thus,  the  3040  cm~^  band  is  the  band  that 
extends  from  3030  to  3050  cm~^. 

Figure  1  shows  our  initial  comparison  of  line  by  line  and  band  model 

transmittances  in  the  Interval  from  2800  to  3400  cm~^.  The  optical  path 

length  in  this  calculation  is  30  km.  A  convenient  measure  of  the  band  model 

error  is  the  logarithm  of  the  ratio  of  the  band  model  transmittance  T  to  the 

_  *  ® 
mean  line  by  line  transmittance  T 

L> 

E  -  LogjQ(Tg/Tj^)  (10) 

This  quantity  is  plotted  in  Fig.  lb.  For  transmittances  generally  in  the 
range  T  >  0.1,  we  see  that  the  band  model  does  quite  well.  However,  when 
T  <  0.1  the  band  model  can  be  in  error  by  several  orders  of  magnitude.  In 
particular,  the  errors  in  the  bands  centered  at  3040  and  3060  cm  ^  are  large 
and  opposite  in  direction.  These  two  are  taken  as  representative  of  low 
transmittance  bands  and  will  be  examined  in  detail  in  the  remainder  of  this 
report. 

The  band  model  and  line  by  line  transmittances  in  the  3040  and  3060  cm~^ 
bands  are  plotted  as  a  function  of  distance  in  Figs.  2a  and  3a  with  the 
corresponding  error  curves  plotted  in  Figs.  2b  and  3b.  In  the  3040  cm~^  band, 
the  error  Is  less  than  unity,  decreases  to  a  minimum,  and  then  begins  to 
Increase,  whereas  for  the  3060  cm~^  band  the  error  Increases  monotonies lly. 


*In  this  report  the  unit  of  error  is  called  an  "order  of  magnitude."  Thus, 
E  “  2  is  a  two  order-of-magnltude  error. 
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These  two  examples  are  representative  of  the  two  basic  shapes  for  error 
curves.  The  reasons  for  these  two  basic  shapes  will  be  discussed  In  more 
detail  In  Section  III-C. 

B.  HoO  LINE  PARAMETER  DISTRIBUTIONS 


The  first  question  that  must  be  answered  concerns  the  line  parameter 
distributions.  Are  the  line  parameters  actually  distributed  In  reasonable 
agreement  with  the  assumptions  of  statistical  band  theory?  In  order  to  answer 
this,  we  have  plotted  the  distributions. 

Figure  A  Is  a  schematic  plot  of  the  lines  In  the  two  bands.  The  height 
of  each  line  Is  the  dimensionless  strength  parameter  ■  S^/S  where  S  Is 
the  average  strength  for  the  band.  In  statistical  band  theory  the  lines  are 
assumed  to  be  randomly  positioned  In  the  band,  and  the  line  strengths  are 
assumed  to  be  distributed  with  a  probability  density  given  by  the  exponential- 
tailed  Inverse  function.^  This  distribution  function  Is 


P(s)  - 


S  In  R 


[e,p  (^)  -  .Kp  (|i)] 


Where  R  and  are  parameters.  It  Is  convenient  to  work  with  the  dimensionless 
strength  parameter  a  ■  S/S  (where  S  Is  the  average  of  the  distribution).  The 
probability  distribution  for  a  Is  easily  shown  to  be 

G(a)  -  S  P(S  o)  (12) 

which  evaluates  to 


G(a)  ~  g  Yn  (^  (“Ao)  -  exp  (-RAo)] 


where 


A  -  S/S 


.  R  ~  1 

m  R  In  R 


This  function  is  represented  by  the  dashed  curves  in  Fig.  5.  In  both 
graphs,  R  ■  10^.  The  actual  density  distributions  of  S/S  in  the  bands  are 
represented  by  the  histograms  in  Fig.  5.  The  function  and  histograms,  in  each 
case,  appear  to  be  in  reasonable  agreement,  vindicating  the  use  of  the 
exponential-tailed  Inverse  distribution. 

The  line  positions  are  assumed  to  be  random.  If  this  is  the  case,  then 
the  spacing  between  lines  has  a  probability  density  given  by 

P(6)  -  i  exp  (-6/6)  (15) 

5 

where  6  is  the  average  value  for  the  spacing  given  by 


The  comparisons  of  this  theoretical  distribution  with  the  actual  histograms 
are  shown  in  Fig.  6.  Again  the  agreement  seems  reasonable. 

In  the  band  model,  the  line  width  is  assumed  to  be  a  constant  equal  to 
the  average  value  6.  The  actual  distributions  of  widths  in  the  3040  and  3060 
cm~^  bands  are  plotted  as  histograms  in  Fig.  7.  The  dashed  curve  is  a 
Gaussian  distribution  function  with  the  same  average  and  variance  as  the 
actual  width  distribution.  It  is  obvious  that  the  widths  are  not  distributed 
normally. 

Finally,  the  absorption  coefficient  k(v),  defined  by  Eq.  (2),  is  plotted 
for  our  two  representative  bands  in  Fig.  8. 

C.  COMPUTER  GENERATED  LINE  PARAMETERS 

Additional  tests  of  the  band  model  were  done  using  a  set  of  computer 
generated  line  parameters  Instead  of  the  experimental  H2O  parameters  obtained 
from  the  line  atlas. 

The  computer  generated  line  parameters  are  random  samples  drawn  from 
infinite  parent  populations  which  are  defined  by  their  probability  distribu¬ 
tion  functions.  The  line  strength  population  is  defined  by  the  exponential- 


15 


Line  Strength  Distribution.  Dashed  curve  is  the  exponential 
tailed  Inverse  function. 


Called  Inverse  distribution,  the  line  wave  number  population  Is  defined  to 
have  a  uniform  distribution,  and  the  line  widths  are  all  set  equal  to  a  con¬ 
stant. 

By  using  these  computer  generated  line  parameters  we  refine  our  analysis 
of  the  band  model  error  since  we  are  now  assured,  as  much  as  possible,  that 
the  line  parameter  distributions  are  In  agreement  with  the  theoretical  assump¬ 
tions.  In  addition,  we  can  study  Che  effect  of  random  fluctuations  on  the 
error  by  generating  many  sets  of  line  parameters  and  calculating  the  error 
curve  for  each. 

A  spectrum  of  250  lines  in  the  wave  number  range  0  to  50  cm  ^  were 
generated.  The  mean  CransmlCtances  were  computed  for  the  20  cm”^  wide  band 
extending  from  15  to  35  cm  The  wave  numbers  for  the  lines  were  computed 
using  the  simple  formula 


*  50  .  (17) 

where  the  are  random  numbers  distributed  uniformly  in  the  interval  0  to 
1.  This  array  of  wave  numbers  was  then  rearranged  so  that  they  were  in  mono- 
tonically  Increasing  order  with  respect  to  the  index  1. 

The  line  width  was  set  equal  to  a  constant  value, 

-  0.075  (18) 

This  constant  Is  approximately  the  same  as  the  average  line  width  In  the  two 

H2O  bands  considered  previously  (see  Fig.  7). 

The  procedure  for  generating  line  strengths  is  more  involved.  The 
relation  between  the  random  variable  X  which  is  uniformly  distributed  in  the 
range  0  to  1,  and  the  relative  line  strength  distribution  given  by  Eq.  (13)  is 

X  -  /  G(o')  da'  -  H(a)  (19) 

0 
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Then  inverting  this  equation  and  using  ■»  we  obtain 

-  S  h"^  (X^)  (20) 

where  the  are  random  numbers,  which  are  independent  of  the  random  numbers 
used  to  generate  the  wave  numbers,  and  is  the  inverse  function  of  H.  The 
functions  H  and  must  be  computed  numerically.  The  line  strengths  computed 
this  way  will  have  an  exponential-tailed  Inverse  distribution  with  an  average 
value  (Although  ?  is  the  average  value  for  the  parent  population,  the 

average  line  strength  of  the  finite  random  sample  that  we  will  compute  will 
deviate  from  this  value.) 

A  different  set  of  line  parameters  will  be  generated  each  time  this 
process  is  carried  out,  since  different  random  numbers  are  used  in  each  run. 
Repeating  the  calculation  many  times  will  generate  an  ensemble  of  line  para¬ 
meter  sets.  (Each  set  is  a  random  sample  drawn  from  the  same  parent  popula¬ 
tion.)  Line  by  line  and  band  model  calculations  were  carried  out  using  these 
computer  generated  line  parameters.  The  results  of  a  typical  calculation  are 
shown  in  Figs.  9  through  13.  These  calculations  were  then  repeated  20  times, 
using  different  random  sample  line  parameters  in  each  calculation.  The  20 
error  curves  are  all  plotted  on  the  same  graph  in  Fig.  14. 

The  error  seems  to  be  composed  of  two  parts,  a  random  fluctuating 
component  and  systematic  component.  The  random  component  dominates  at  inter¬ 
mediate  distances  and  shows  no  bias,  i.e.,  it  is  just  as  likely  positive  or 
negative.  The  systematic  component  dominates  at  large  distances  and  increases 
without  limit  as  the  distance  Increases.  This  just  means  that,  in  this  limit, 
the  band  model  always  predicts  larger  transmlttances  than  the  line  by  line 
calculation.  Such  behavior  is  easily  understood  by  examining  the  asymptotic 
behavior  of  the  band  model  and  line  by  line  transmlttances.  The  asymptotic 
form  for  the  band  transmittance,  obtained  from  Eq.  (5),  is 

T(x)  exp[-2  (^  x)^^^]  (21) 

e 
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The  line  by  line  band  transmittance  Is  given  by  Eq.  (4).  It  Is  obvious 
that  when  x  is  very  large  only  the  values  of  k(v)  in  the  vicinity  of  the 
minimum  value  will  contribute  to  the  integral.  The  Integral  can  be 

evaluated  by  the  steepest  descent  method  to  obtain  the  asymptotic  formula 

T  -*•  (-rS— )  exp(-k  x)  (22) 

k  X  min 

where  k"  is  the  second  derivative  of  k(v)  evaluated  at  the  minimum.  The 
asymptotic  form  of  the  error  E  is  then  derived  from  Eq.  (10), 

E  .  0.434  x]  (23) 

The  essentially  linear  Increase  at  large  distances  is  evident  in  all  the  error 
curves  plotted  in  Fig.  14. 

The  random  component  of  the  error  will  be  discussed  in  Section  III-C. 


III.  NEW  TRANSMITTANCE  APPROXIMATIONS 


The  basic  result  of  statistical  band  model  theory  Is  the  formula,  Eq. 
(5),  for  computing  the  mean  band  transmittance.  This  formula  has  three 
adjustable  band  model  parameters,  which  are  determined  either  by  fitting  to 
experimental  data  or  by  calculating  directly  from  line  parameter  data  using 
Eqs.  (6)  through  (8).  The  failure  of  this  method  for  moderate  to  long  optical 
path  distances  has  been  demonstrated  In  the  previous  section.  In  this  section 
we  derive  alternative  parameterized  formulae  for  computing  T(x)  which  are 
accurate  for  all  values  of  the  optical  path  length. 

Rewrite  the  Integral  In  Eq.  (4)  as  a  simple  numerical  quadrature 

1  ” 

T(x)  *  —  E  expI-k(Vj^)x]  6v  (24) 

where  N  Is  the  number  of  quadrature  points,  6v  Is  the  spacing  between  points, 
and  Av  Is  the  bandwidth.  Since  Av  ■  N6v,  this  can  also  be  written 

1  N 

“  M  ^  exp[-k(v  )  x]  (25) 

N  1.1  1 

N  must  be  large  enough  to  ensure  adequate  accuracy.  (For  the  calculations  In 
the  previous  section  we  used  N“1000.)  The  numerical  ordering  of  the  terms  In 
the  sum  does  not  matter.*  Thus  the  array  k(Vj^)  of  discrete  k  values  can  be 
rearranged  In  monotonlcally  Increasing  order.  The  points  were  originally 
spaced  6v  “  Av/N  units  apart.  The  rearranged  points  are  spaced  ^  *  1/N  units 
apart  In  the  unit  Interval  0  ^  p  ^  1.  It  Is  useful  to  regard  these  points  as 
defining  a  monotonlcally  Increasing  continuous  function  of  p  in  this 
Interval.  (One  could  define  this  function,  for  example,  by  connecting  the 


The  basic  Idea  of  reordering  k  values  is  quite  old.  Application  of  the 
method  and  references  to  Its  previous  use  are  given  In  Ref.  5. 

^A.  Arklng  and  K.  Grossman,  J.  Atmos.  Scl.  29,  937  (1972). 
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points  with  straight  line  segments.)  We  call  this  function  k(p)  the  "mono- 
tonic  absorption  function."  The  formula  for  the  transmittance  can  be 
expressed  In  terms  of  k(p) 

_  1 

T(x)  ■  /  exp[-k(p)  x]  dp  (26) 

0 

or.  In  discrete  form 

1  ” 

T(x)  - I  exp(-k(p. )x]  (27) 

1-1 

The  function  k(p)  Is  shown  for  the  3040  and  3060  cm~^  bands  of  H2O  In 
Fig.  15.  These  should  be  compared  to  Fig.  8  where  the  absorption  coefficients 
are  shown  plotted  as  k(v)  In  their  natural  order.  It  Is  obvious  that  the 
k(v)  functions  cannot  be  approximated  by  any  simple  analytic  function.  The 
monotonlc  functions  k(p)  however  may  be  amenable  to  simple  analytic  approxima¬ 
tions.  Both  of  the  graphs  of  k(p),  shown  In  semi-log  plots  In  Fig.  15,  appear 
to  be  roughly  linear.  Therefore  the  first  approximation  we  will  try  Is  just  a 
simple  exponential  function 

k(p)  =  kg  exp(bp)  (28) 

A.  TWO-PARAMETER  APPROXIMATION 

In  the  previous  section  It  was  shown  that  the  long  range  behavior 
of  T(x)  Is  dominated  by  k^^ Since  we  want  the  long  range  behavior  to  be 
correct,  we  define  Icq  -  k^^  ^  «  k(0).  At  the  other  extreme,  the  very  short 
range  behavior  of  T(x)  Is  determined  by  the  average  value  of  k.  This  Is 
easily  proven.  For  very  small  values  of  x,  the  exponential  In  Eq.  (26)  can  be 
replaced  by  the  first  two  terms  of  Its  power  series  expansion,  thus 
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Fig.  15.  Monotonlc  Absorption  Function  k(p).  Comparison  of  line  by 
line  and  two-parameter  model.  Tick  marks  are  upper  Integra 
tlon  limits  used  in  Eq.  (26)  for  90%  accuracy  at  distances 
Indicated  In  kilometers. 


<k>  “  J  k(p)  dp 
0 


is  the  average  value  of  k. 


The  parameter  b  in  Eq.  (28)  is  chosen  so  that  the  average  value  of  the 
exponential  approximation  is  equal  to  <k>.  Thus ,  we  obtain 


<k>  =  k  y  exp(bp)  dp 
^  0 


b  1 

e  -  1 


which  can  be  solved  numerically  for  b.  Using  these  values  for  the  parameters 
kQ  and  b,  the  exponential  approximation  to  k(p)  is  plotted  as  the  dashed 
straight  lines  in  Fig.  15.  The  approximate  transmittance  is  then  computed  by 
substituting  the  function  k(p)  given  by  Eq.  (28)  into  Eq.  (26)  and  integrating 
the  resulting  expression.  The  integral  can  be  evaluated  analytically. 

To  accomplish  this,  the  variable  of  integration  in  Eq.  (26)  is  changed 
from  p  to  k 


T(x)  ■  f.  explkx]  f(k)  dk 


where  the  function  f(k)  is 


dp(k) 

dk 


kj  «  kp  exp(b) 


[The  meaning  of  the  function  f(k)  will  be  discussed  in  more  detail  later.] 


Equation  (28)  Is  easily  Inverted  to  obtain 

p(k)  - -g  InCk/kp)  (36) 

and  therefore,  using  Eq.  (3A) 

£(k)  -  -L  (37) 

Substituting  this  expression  Into  Eq.  (33)  gives 

T(x)  -  i  [E^  (k^  x)  -  E^  (k^  x)J  (38) 

where  E2(x)  Is  the  exponential-integral  function  defined  by^ 

«0  “t 

Ej(x)  -  ;  dt  (39) 

X 

Very  efficient  methods  are  available  for  the  numerical  evaluation  of  the 
exponential  Integral  function.^  The  function,  Eq.  (38),  has  three  parameters 
kQ,  kj,  and  b.  However  only  two  are  Independent.  Using  Eq.  (35)  we  express  b 
In  terms  of  kQ  and  kj^  obtaining  our  final  two-parameter  expression  for  T(x) 

'‘O’ 

This  function  was  evaluated  numerically  and  the  approximate  mean 
transmittance  curves  are  shown  plotted  In  Figs.  16a  and  17a  along  with  the 
precise  line  by  line  results  for  comparison.  Below  these  graphs  In  Figs.  16b 
and  17b  are  the  error  curves  for  the  two-parameter  approximation  and  also,  for 
comparison,  the  error  curves  for  the  band  model  transmlttances  (see  Figs.  2b 
and  3b). 


^M.  Abramowltz  and  I.  Stegun,  Handbook  of  Mathematical  Functions 
Dover,  New  York  (1965) 


Fig.  17.  Comparison  of  Line  by  Line  (T  )  and  Two-Parameter_Mode^l  (T-) 
Transmlctances  for  3060  cin“^  Band.  Error  ■  Log^gd^/T^).  Th 
band  model  error  is  also  shown  for  comparison. 


The  range  of  integration  in  Eq.  (26)  for  the  variable  of  integration  p  is 
0  ^  p  ^  1.  For  large  values  of  the  distance  x,  the  contribution  from  the 
upper  portion  of  this  integration  range  is  very  small.  In  order  to  demon¬ 
strate  this  quantitatively,  we  have  computed  the  upper  limit  of  integration 
required  to  compute  90%  of  the  value  of  T(x)  for  various  fixed  values  of  x. 
These  upper  integration  limits  are  shown  as  tick  marks  in  Fig.  15  and  are 
labeled  with  the  appropriate  value  of  distance  x  in  kilometers.  Ue  see  that 
even  for  fairly  small  distances,  the  upper  portion  of  the  integration  range 
does  not  make  much  contribution.  This  would  suggest  that  any  analytic  fit  to 
k(p)  should  be  weighted  to  have  the  least  error  in  the  lower  part  of  the  p 
range.  The  present  procedure  for  calculating  the  slope  parameter  b  overempha¬ 
sizes  the  large  k  values.  From  Fig.  IS  we  see  that  a  better  fit  would  be 
obtained  in  the  lower  portion  of  the  k(p)  curve  if  the  slope  parameter  b  were 
less. 

A  new  procedure ,  which  gives  more  weight  to  the  lower  k  values ,  was  tried 
for  fitting  the  analytic  function,  Eq.  (28),  to  k(p).  The  parameter  kQ  is 
still  defined  to  be  the  minimum,  kQ  *  k(0).  However,  instead  of  computing  the 

average  value  of  k  as  in  Eq.  (30),  we  now  compute  the  average  value  of  the 

natural  log  of  k 

1 

<  £n(k)  >  *  /  £n  lk(p)]  dp  (41) 

0 

The  parameter  b  is  determined  by  requiring  that  the  average  value  of  the 

natural  log  of  the  exponential  approximation,  Eq.  (28),  is  equal  to  <  £n(k)  > 

1 

<  jin(k)  >  =  /  £n  [k  exp(bp)]  dp  (42) 

0  ^ 

This  Is  easy  to  solve  and  we  obtain 

b  =  2  l<£n(k)>  -  £n(kQ)J  (43) 

The  new  value  for  the  parameter  k^  is  then  computed  by  substituting  this  value 
of  b  into  Eq.  (35)  and  the  transmittance  is  computed  using  Eq.  (40). 


This  new  analytic  approximation  to  k(p)  is  shown  plotted  as  the  dashed 
straight  line  in  Fig.  18.  The  slope  of  the  line  has  been  reduced  from  that  of 
Fig.  15  and  a  better  fit  to  k(p)  is  obtained  in  the  lower  range  of  the 
curve.  In  Figs.  19  and  20,  which  show  the  transmittances  and  error  curves  for 
our  two  example  bands,  the  error  has  been  reduced  from  that  shown  in  Figs.  16 
and  17.  Thus,  the  second  method  for  computing  b,  which  gives  more  weight  to 
the  lower  k  values,  is  slightly  superior.  It  has  a  very  simple  graphical 
interpretation  in  the  semi-log  plots  shown  in  Fig.  18.  The  area  under  the 
straight  line  approximation  is  equal  to  the  area  under  the  line  by  line 
k(p).  (Unless  otherwise  stated,  in  any  future  reference  to  the  two-parameter 
method,  the  parameters  are  computed  by  the  <  £n(k)  >  method.) 


B.  THREE-PARAMETER  APPROXIMATION 


Any  further  significant  improvement  in  accuracy  can  only  be  accomplished 
by  increasing  the  flexibility  of  the  analytic  function  used  to  approximate 
k(p).  We  have  done  this  by  dividing  the  integration  range  into  two  parts,  0  _< 
p  ^  1/2  and  1/2  ^  1,  and  approximating  k(p)  in  each  of  these  regions  by  an 
exponential  function 


k(p) 


Rq  exp[bjpJ ,  0  <  p  <  1/2 

ki/2  explb2(p-l/2) j ,  1/2  <  p  <  1 


(44) 


The  function  is  required  to  be  continuous  which  implies 

ki/2  =  k^  exp[b^/2]  (45) 

The  maximum  value  of  this  function  is 


The  inverse  function  to  Eq 


p(k) 


ki  -  kj^2  explb2/2] 

(44)  is 

£n(k/kQ)  kp  <  k  <  k^/2 
£n(k/kj)  k^^2  <  k  <  k^ 


(46) 


(47) 


and  therefore  from  Eq.  (34) 


3040  cm  *  BAND 


Icq  <  k  <  kj^2 


(48) 


f(k) 


kj/2  <  k  <  kj 


Note  that  f(k)  Is  discontinuous  at  k  *  kj^2* 

The  minimum  value  of  the  analytic  approximation  is  kQ.  Just  as  in  the 
previous  cases,  this  is  defined  to  be  equal  to  the  true  minimum  value  k^^^  of 
k(p), 


min 


k(0) 


(49) 


The  parameters  and  b2  are  defined  similarly  to  the  previous  case.  If  we 
calculate  the  average  values  of  lnlk(p)J  In  each  half  region 

1/2 

<  Jln(k)  >1  *  /  Jin[k(p)l  dp  (50) 

0 

1 

<  in(k)  >,  *  /  ln(k(p)]  dp  (51) 

1/2 

and  equate  these  to  the  values  obtained  when  the  analytic  approximation,  Eq. 
(44),  Is  substituted  for  k(p)  In  Eqs.  (50)  and  (51), 

b^  «  4  [2  <ln(k)  >j  -  indc^)]  (52) 

and 

b^  -  4(2  <  in(k)  >2  -  in(kj^2^J 

The  value  of  kj^2  ^9*  (53)  is  computed  using  Eq.  (45).  With  all  the  para¬ 
meters  defined,  the  analytic  approximation,  Eq.  (44),  Is  evaluated  and  plotted 
in  Fig.  21.  This  should  be  compared  with  the  previous  two-parameter  approxi¬ 
mation  In  Fig.  18. 
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The  transmittance  is  obtained  by  substituting  the  expression  for  f(k) 
given  by  Eq.  (A8)  into  Eq.  (33)  and  evaluating  the  integrals.  The  result  is 


T(x) 


^7 


x)  - 


Ei(ki/2  x)l 


x)  - 


Ei(ki 


X)] 


(54) 


Using  Eqs.  (45)  and  (46)  we  can  express  bj  and  b2  in  terms  of  kQ,  kj^y2* 
kj.  The  final  form  for  our  three-parameter  formula  is 


T(x) 


_ 1 _ 

2  £n(k^^2/^0^ 


(E^(kQ  x)  -  E^(kj^2 


2  An(kj^/kjy2 


(55) 


Transmittances  computed  using  this  approximation  are  shown  plotted  in 
Figs.  22  and  23  along  with  the  error  curves.  As  can  be  seen,  this  three- 
parameter  approximation  is  excellent  over  the  entire  range  of  distances. 

In  the  discussion  presented  so  far,  the  two-parameter  and  three-parameter 
approximations  have  been  presented  as  approximations  to  line  by  line  calcula¬ 
tions.  Obviously,  they  can  also  be  used  as  a  convenient  fit  to  experimental 
transmittance  data. 

C.  k-DISTRIBUTION  FUNCTION 

The  function  f(k)  is  defined  by  Eq.  (34).  In  differential  form  it  is 


dp  ■  f(k)  dk 


(56) 


Integrating  this  expression  gives 

kp 

Ap  =  P2  -  Pj  “  f(k)  dk  (57) 

The  quantity  Ap  is  the  fraction  of  the  spectral  interval  for  which  k  is  in  the 
range  ^  k  _<  k2.  If  k  is  considered  to  be  a  random  variable,  then  f(k)  is 
the  probability  density  distribution  of  k  and  Ap  is  the  probability  that  a 
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randomly  chosen  k  will  be  in  the  interval  k  _<  k2.  For  this  reason  the 

function  f(k)  has  been  called  by  Arking  and  Grossman^  the  k-distribution 
function. 

Given  the  k-distribution  function  f(k),  the  mean  transmittance  is  given 
by  Eq.  (33).  In  this  expression  the  lower  and  upper  integration  limits  have 
been  expressed  as  the  finite  values  kQ  and  kj^.  This  is  in  conformity  with  our 
expectations  for  real  spectra  where  kQ  >  0  and  k^  <  «.  It  will  be  useful  now 
to  extend  these  values  to  their  ultimate  limits  and  write 

T(x)  ■  j  exp[k  x]  f(k)  dk  (58) 

0 

[This  more  general  formula  is  true  even  for  real  spectra,  since  we  only  need 
to  define  f(k)  ■  0  outside  the  range 

Domoto^  recognized  that  this  relation  defines  T(x)  as  the  Laplace 
transform  of  f(k)  and,  conversely,  f(k)  is  the  inverse  Laplace  transform  of 
T(x).  He  applied  this  procedure  to  the  statistical  band  model  where  T(x) 

is  given  by  Eq.  (5).  The  Inverse  transform  he  obtained  is 

_  3  1/2 

f<k)-|  .xp[a(2  (59) 

Iv  Iv 

where  a  *  This  distribution  is  defined  on  the  entire  interval  0  ^  k 

K.  *0. 

The  k-distributions  f(k)  for  the  two  and  three  parameter  models  are  given 
by  Eqs.  (37)  and  (48),  respectively.  The  exact  k-distrlbution  can  be  computed 
numerically  from  a  line  by  line  calculation.  These  functions  are  plotted  in 
Figs.  24  and  25.  Figure  24  compares  a  band  model,  line  by  line  and  the  two- 
parameter  model.  Figure  25  is  the  same,  but  plots  f(k)  for  the  three- 
parameter  model  instead  of  the  two-parameter  model.  Note  that  the  three- 
parameter  f(k)  is  discontinuous. 


^G.  Domoto,  J.  Quant.  Spectrosc.  Radiat.  Transfer  14,  935  (1974). 


25.  k-Dlstrlbutlon  Functions.  Same  as  Fig.  24  except  that  the 
solid  line  Is  the  three-parameter  model. 


The  cumulative  distribution  function  is  obtained  by  integrating  the 


probability  density  function 


f(k’)  dk*  -  F(k) 


The  inverse  of  this  cumulative  probability  function  is  the  function  we  have 
previously  defined  as  the  "monotonic  absorption  function"  k(p).  Thus 

k(p)  -  f"^(p)  (61) 

This  procedure  provides  a  means  of  computing  the  monotonic  absorption  function 
k(p)  for  the  statistical  band  model.  Both  steps  of  the  procedure  must  be  done 
numerically:  first  a  numerical  Integration  of  f(k)  to  obtain  F(k)  and  then  a 
numerical  interpolation  to  obtain  F  ^(p).  The  results  are  plotted  in  Fig.  26. 
The  exact  results  are  also  plotted  for  comparlsdn.  This  graph  should  be  com¬ 
pared  with  Figs.  15,  18,  and  21,  where  k(p)  for  the  two-  and  three-parameter 
models  are  plotted. 

The  functions  T'(x)  and  k(p)  are  transforms  of  each  other.  One  can  be 

computed  from  the  other.  The  function  k(p)  is  calculated  from  T(x)  by  the 

procedure  just  outlined  and  T(x)  is  calculated  from  k(p)  by  Eq.  (26).  The 
behavior  of  T(x)  for  large  values  of  x  is  determined,  for  the  most  part,  by 
the  values  of  k(p)  in  the  interval  near  p  0.  This  is  Illustrated  in  Fig.  15 
where  the  upper  limit  of  Integration  used  in  Eq.  (26)  to  calculate  90%  of  the 
final  value  of  T(x)  is  plotted  for  various  values  of  x.  In  the  limit  as 
X  -»•  CO,  T(x)  is  determined  by  the  single  point  at  p  ■  0,  k(0)  “  ^in 
(22) J.  Thus,  in  order  to  compute  T(x)  accurately  for  large  optical  path 

lengths,  it  is  necessary  to  have  an  accurate  approximation  of  the  function 

k(p)  at  and  near  p  ■  0. 


’S.  L.  Meyer,  Data  Analysis  for  Scientists  and  Engineers,  Wiley,  New  York 
(1975)  p.  20,  103. 
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The  two-  and  three-parameter  approximations  we  have  proposed  take  advan¬ 
tage  of  the  relatively  simple  form  of  the  function  k(p)  by  approximating 
ln(k(p)]  by  one  or  two  straight  line  segments.  These  approximations  are 
constrained  to  be  equal  to  '^mln  at  p  *  0.  This  ensures  that  our  approximation 
will  have  the  correct  long  range  behavior. 

The  statistical  band  model,  on  the  other  hand,  is  basically  limited  In 

accuracy  for  large  values  of  x  because  the  Information  contained  In  the  band 

model  parameters  k,  y ,  and  6^  Is  not  sufficient  to  determine  More 

generally  It  Is  not  sufficient  to  determine  k(p)  in  the  interval  near  p  =  0. 

The  best  that  can  be  done  with  this  Information  is  to  determine  a  distribution 

of  k  (p)  functions  which  are  compatible  with  the  parameters.  (The  subscript 
a 

a  labels  the  functions  In  this  distribution.)  The  functions  k  (p)  determine 

_  a 

a  family  of  mean  transmittance  functions  T  (x).  These  in  turn  determine  a 

a 

family  of  error  curves  E  (x),  which  will  look  very  much  like  the  family  of 
error  curves  shown  In  Fig.  14. 

If  the  only  Information  we  are  given  about  a  band  is  the  statistical  band 
model  parameters,  then  each  of  the  error  curves  equally  probable. 

The  error  curve  that  applies  In  any  particular  case  can  be  thought  of  as  just 
a  random  selection  from  this  family  of  curves. 

Before  leaving  this  section.  It  Is  Interesting  to  examine  the  functions 
k(p)  plotted  In  Fig.  26  In  more  detail.  For  convenience,  the  band  model  func¬ 
tion  will  be  written  kg(p)  and  the  line  by  line  function  will  be  written  as 
kj^(p).  The  general  behavior  of  ^^g(p)  is  that  It  turns  sharply  downward  near 
p  “  0  and  approaches  the  value  kg(0)  “  0,  which  Is  always  less  than  = 

kj^(O).  This  Is  a  manifestation  of  the  Incorrect  asymptotic  form  of  the  band 
model  transmittance  function  given  by  Eq.  (21). 


The  family  of  error  curves  shown  in  Fig.  15  Is  not  exactly  the  same  as  those 
described  here.  In  Fig.  15  the  parent  populations  were  fixed,  whereas  In  the 
case  decrlbed  here  the  band  model  parameters  k,  Y,  and  6g  are  held  constant. 
After  the  calculations  for  this  report  were  complete,  we  discovered  a  simple 
method  to  generate  line  parameters  that  have  fixed  band  model  parameters. 
However,  we  believe  that  the  family  of  error  curves  shown  In  Fig.  15  Is  at 
least  a  good  qualitative  and  also  a  semiquantltatlve  picture  of  the  behavior 
that  would  result  when  the  band  model  parameters  are  held  fixed. 


49 


In  the  region  0  _<  p  _<  P2  where  P2  ~  0.5,  the  relationship  between  kg(p) 

and  kj^(p)  will  generally  be  one  of  two  types.  Either  kg(p)  <  kj^(p)  as  in  the 

lower  graph  of  Fig.  26  or  k^Cp)  and  kL(p)  will  Intersect  at  some  small  value 

Pj  with  kg(p)  >  kj^(p)  in  the  region  Pj  ^  P  ^  P2  as  in  the  upper  graph  of  Fig. 

26.  For  the  first  case,  we  can  predict  from  Eq.  (26)  that  T„(x)  >  T,  (x)  and 

B  L 

thus  the  error  ratio  will  be  greater  than  1.  This  behavior  is  Illustrated  In 
Fig.  3b. 

In  the  second  case,  the  behavior  is  more  complicated.  For  moderate 
values  of  x  the  result  will  be  T  (x)  <  T-(x).  However,  as  x  Increases,  the 
effective  range  of  integration  in  Eq.  (26)  decreases.  At  some  point  it  will 
fall  entirely  within  the  Interval  0  P  ^  Pj  where  kg(p)  <  kL(p).  When  this 
occurs  we  obtain  T  (x)  >  T  (x).  Thus  the  functions  T,  (x)  and  T„(x)  eventually 
cross  and  the  error  ratio  varies  from  less  than  1  to  greater  than  1.  The 
beginning  of  this  behavior  is  Illustrated  in  Fig.  2b  which  presumably  would 
follow  the  scenario  just  outlined  if  x  were  extended  beyond  100  km.  The  error 
curves  shown  in  Fig.  14  Illustrate  both  types  of  behavior. 

D.  MONTE  CARLO  METHOD 

The  term  "Monte  Carlo”  usually  refers  to  a  computational  procedure  in 
which  a  large  (or  infinite)  distribution  of  values  of  some  quantity  or 
parameter  is  replaced  by  a  manageably  small,  unbiased  random  sample  of  the 
distribution.  For  example,  a  physical  quantity  may  be  the  average  value  of 
some  distribution.  The  Monte  Carlo  approximation  is  computed  by  generating  a 
random  sample  of  the  distribution  and  then  averaging  the  random  sample.  In 
the  present  case,  the  Monte  Carlo  method  is  a  practical  procedure  for  calcu¬ 
lating  the  parameters  for  the  two-  and  three-parameter  approximations  from  an 
unbiased  random  sample  of  the  k-dlstrlbutlon. 

The  k-dlstrlbutlon  is  the  continuously  Infinite  collection  (or 
population)  of  k  values  that  are  defined  by  the  function  k(v)  in  the  band 
Interval  Av  without  regard  to  their  order.  An  unbiased  random  sample  of  this 
infinite  population  of  k  values  can  be  generated  by  randomly  selecting  M  wave 


numbers  (i  “  1,  M)  from  the  band  interval  and  then  calculating  k(v^)  using 
Eq.  (2).  The  random  wave  numbers  are  computed  from  random  numbers  in  the 
interval  0  to  1  by  the  formula 

Vf  “  (vy  ~  X^,  i  *  1,  M  (62) 

where  v„  and  v,  are  the  upper  and  lower  boundaries  of  the  band. 

U  Lj 

The  least  value  of  k  in  this  random  sample  is  a  good  approximation  to 

kmin  and  is  set  equal  to  the  parameter  kg  for  both  the  two-  and  three- 

parameter  approximations. 

The  average  value  of  in(k)  for  the  M  values  of  k  in  the  random  sample  is 
the  Monte  Carlo  approximation  to  the  quantity  <jin(k)>  defined  by  Eq.  (41). 
The  parameter  k^  for  the  two-parameter  approximation  is  then  calculated  by 
combining  Eqs.  (35)  and  (43)  to  obtain 

k,  *  ^  exp[2  <tn(k)>]  (63) 

For  the  three-parameter  model  we  must  compute  Monte  Carlo  approximations 
for  the  quantities  <in(k)>^  and  <Jin(k)>2  defined  by  Eqs.  (50)  and  (51), 

respectively.  The  quantity  k(p*l/2)  is  the  median  value  of  k  which  we 

designate  1^^^.  This  is  the  value  of  k  such  that  half  the  elements  of  the 

distribution  exceed  it  in  value  and  half  are  less  in  value.  The 

quantity  <An(k)>^  is  the  average  value  on  tn(k)  for  k  <  and  <£n(k)>2  is 

the  average  value  of  Jin(k)  for  k  >  Thus  in  the  Monte  Carlo 

approximation  we  divide  the  random  sample  of  M  elements  into  two  groups,  each 
with  M/2  elements  such  that  any  k  value  in  the  first  group  is  less  than  any  k 
value  in  the  second  group.  (In  order  to  avoid  any  ambiguities  in  this 

procedure,  we  always  choose  M  to  be  an  even  integer.)  The  average  value 

of  £n(k)  in  the  first  group  is  an  approximation  to  <in(k)>^  and  the  average 
value  of  jin(k)  in  the  second  group  approximates  <Jin(k)>2»  The  parameters  ^2/2 
and  kj  for  the  three-parameter  model  are  then  computed  by  combining  Eqs.  (45), 
(46),  (52),  and  (53)  to  obtain 

k./o  ■ expl4<4n(k)>, ]  (64) 


and 


k  - -j— ^  exp(4<jin(k)>,] 


(65) 


In  Figs.  27  and  28  we  show  the  results  obtained  by  this  method  using  the 
three-parameter  approximation  and  a  random  sample  of  M  ■  50  values  of  k. 
These  results  are  almost  as  good  as  the  results  obtained  previously  using  1000 
points  but  required  about  1/20  the  computational  effort.  The  Monte  Carlo 
procedure  Is  thus  an  efficient  method  for  carrying  out  (approximate)  line  by 
line  calculations. 
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Pig.  27.  Comparison  of  Line  by  Line  (T  )  and  50  Point  Monte  Carlo  (T_p.) 
Transmittances  for  3040  cm“^  Band.  Error  ■  Log,^(f-^/f, ). 
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.Ine  by  Line  (Tj^)  and  50  Point  Monte 
insmlttances  for  3060  cm~^  Band.  Error  “ 


IV.  SUMMARY  AND  DISCUSSION 


In  this  report,  we  have  investigated  the  reliability  of  the  statistical 
band  model  by  comparing  the  model  with  precise  line  by  line  calculations  and 
we  have  also  derived  two  new  nonstatistical  band  model  approximations. 

The  first  part  of  the  study,  which  is  an  evaluation  of  the  statistical 
band  model,  is  contained  mainly  in  Section  II  with  some  additional  discussion 
in  III-C.  The  study  was  carried  out  on  a  system  of  H2O  absorbers  in  the  wave 
number  range  2800  to  3400  cm~^.  Tests  were  also  made  using  sets  of  computer 
generated  line  parameters.  The  results  of  these  calculations  were  presented 
graphically. 

Figures  5  and  6  show  that  the  distribution  of  H2O  line  strengths  and  line 
spacings  are  in  reasonable  agreement  with  the  theoretical  assumptions  made  in 
the  statistical  band  model.  The  line  by  line  and  band  model  txansmlttances 
are  compared  in  Figs.  1,  2,  3,  13,  and  14.  It  is  concluded  from  these  (and 
others  not  presented  here)  that  the  statistical  band  model  is  fairly  reliable 
in  the  short  optical  path  length  regime  in  which  T(x)  >  0.1,  but  for  long 
paths  where  T(x)  <  0.1,  the  transmittances  can  be  in  serious  error  by  orders 
of  magnitude.  This  error  has  a  random  component  that  dominates  at  intermedi¬ 
ate  distances  and  a  systematic  component  that  dominates  at  very  long  dis¬ 
tances. 

The  random  error  arises  because  the  information  contained  in  the  band 
model  parameters  is  not  sufficient  to  define  a  unique  k-dlstrlbution,  but 
rather  is  compatible  with  an  entire  ensemble  of  k-dlstrl but ions  from  which  one 
has  been  randomly  selected. 

The  systematic  error  arises  because  of  certain  simplifying  assumptions 
made  in  the  derivation  of  the  band  model  transmittance.  All  values  of  the 
strength  parameter  and  all  values  of  line  spacing  from  0  to  <»  were  allowed. 
As  a  result,  the  k-dlstrlbutlon  extends  from  0  to  “  whereas  any  real  distribu¬ 
tion  has  finite  limits,  >  0  and  k^^^  <  “.  The  unphysical  k  values  in  the 

range  0  to  k^j^u  result  in  an  erroneous  asymptotic  behavior  for  the  band  model 


transmittance  that  Is  consistently  too  large.  This  systematic  component  of 
the  error  Is  easily  computed  from  the  asymptotic  formulas  and  Is  given  In  Eq. 
(23). 

The  second  part  of  this  report  Is  concerned  with  the  derivation  and 
testing  of  two  new  nonstatlstlcal  band  model  approximations.  They  are 
referred  to  simply  as  the  two-parameter  and  three-parameter  approximations. 
The  formulas  for  these  two  approximations  are  given  by  Eqs.  (40)  and  (55). 
They  are  compared  with  the  exact  results  and  with  the  statistical  band  model 
In  Figs.  19,  20,  22,  and  23.  We  conclude  that  the  two-parameter  model  Is 
sometimes  slightly  Inferior  to  the  statistical  band  model  for  short  optical 
paths  but  Is  always  much  superior  for  long  paths.  The  three-parameter  model 
Is  uniformly  excellent  for  all  path  lengths. 

The  present  study  was  a  preliminary.  Investigation  and  was  limited  In  Its 
scope.  A  more  complete  study  should  repeat  most  of  the  calculations  In  this 
report  over  a  broader  range  of  conditions  Including  much  higher  and  lower 
temperatures  and  pressures,  for  other  portions  of  the  H2O  line  system,  and  for 
several  other  molecules,  especially  C02*  Also  other  line  profiles  should  be 
studied. 

Several  other  topics  that  could  be  Included  In  a  new  study  would  be  a 
study  of  the  temperature  and  pressure  dependence  of  the  model  parameters  kg, 
^1/2'  ^1'  ®  generalization  of  our  new  approximations  to  systems  with  non- 

uniform  temperatures  and  pressures;  and  the  development  of  practical  numerical 
techniques  for  fitting  both  the  two-  and  three-parameter  formulas  to  experi¬ 
mental  transmittance  data. 
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